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Abstract 

Apart from the higher limits of isospin and temperature, the properties of atomic nuclei are in¬ 
triguing and less explored at the limits of lowest but finite temperatures. At very low temperatures 
there is a strong interplay between the shell (quantal fluctuations), statistical (thermal fluctua¬ 
tions), and residual pairing effects as evidenced from the studies on giant dipole resonance (GDR). 
In our recent work [Phys. Rev. C 90 , 044308 (2014)], we have outlined some of our results from 
a theoretical approach for such warm nuclei where all these effects are incorporated along within 
the thermal shape fluctuation model (TSFM) extended to include the fluctuations in the pairing 
field. In this article, we present the complete formalism based on the microscopic-macroscopic 
approach for determining the deformation energies and a macroscopic approach which links the 
deformation to GDR observables. We discuss our results for the nuclei 97 Tc, 120 Sn, 179 Au, and 
208 Pb, and corroborate with the experimental data available. The TSFM could explain the data 
successfully at low temperature only with a proper treatment of pairing and its fluctuations. More 
measurements with better precision could yield rich information about several phase transitions 
that can happen in warm nuclei. 
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I. INTRODUCTION 


The study of giant dipole resonance (GDR) at high temperature (T) and angular mo¬ 
mentum (J) has been an interesting area of research which has revealed several structural 
properties of nuclei at extreme conditions. Being a fundamental mode of photo excitation, 
GDR can probe nuclei at extreme conditions and even those with exotic structures [1] . Sev¬ 
eral earlier works on GDR focused on the liigh-J regime m whereas recent studies at 
extreme isospins have potential astrophysical implications |1H3- The GDR at low T is also 
relatively less explored and its studies have gained acceleration in recent times MB- Ex¬ 
perimentally, it is very difficult to populate nuclei at low excitation energies, but it is still 
feasible due to recent developments in the experimental facilities. Several properties of nuclei 
at low T are still not clear: for example, the existence of pairing phase transition, the order 
of it if it exists, the role of fluctuations, etc. The low-T region is quite intriguing because the 
microscopic effects such as shell (quantal) and pairing effects are quite active, although they 
are modified by thermal effects. Since the nucleus is a tiny system, the thermal fluctuations 
inherent in finite systems are expected to be large. The shape degrees of freedom being 
crucial for nuclear structure, the deformation parameters are closely associated with the 
order parameters for the related transitions. Hence the thermal shape fluctuations (TSF) 
are the most dominant ones, and at low T the fluctuations in the pairing field can also 
contribute significantly. Many models [T21 ITflj have been used to study the effect of both 
these fluctuations separately, but the combined effect of these two was not investigated until 
our recent efforts PUCE]. 

Various theoretical approaches have been introduced to investigate the GDR. In a macro¬ 
scopic approach, GDR is a collective mode of excitation of nuclei caused by the out-of-phase 
oscillations between the proton and neutron fluids under the influence of the electromagnetic 
field induced by an emitted/absorbed photon. Here, the GDR couples directly to the shape 
of the nucleus hence providing corresponding structure information. The thermal shape 
fluctuation model (TSFM) is based on such a macroscopic approach and explains the in¬ 
crease in GDR width with T by taking in to account the average of GDR cross sections over 
the shape degrees of freedom. The TSFM is quite successful in interpreting several GDR 
measurements at high T and J and hence is often used by experimentalists as theoretical 
support for their data. 
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In an alternative approach, the GDR observables can be calculated utilizing a linear 
response theory incorporating the TSF within the static path approximation mm- Other 
rigorous microscopic approaches exist [2T, 22] but are mostly limited to the study of low- 
lying states. In microscopic approaches, such as the phonon-damping model (PDM), GDR 
damping is explained through coupling of the GDR phonon to particle-hole, particle-particle, 
and hole-hole excitations [8J m I23H26]. The PDM can explain the increase in GDR width 
with T, at moderate T and the width saturation at high T. 

Apart from the above discussed models, a few phenomenological parametrizations have 
been reported HDU2Z], which are very successful in explaining the global trend of the GDR 
width as a function of T and J. These parametrizations are constructed to mimic the 
results of TSFM and hence have only empirical basis without any microscopic or macroscopic 
foundation. However, in recent literature mum Eli Eg, these parametrizations are referred 
synonymously to models and hence should not be understood as a variant of TSFM. 

The damping of GDR width at low T and the strong influence of pairing on GDR width 
were first suggested by one of the present authors [ 8 ]. The first low-T GDR width mea¬ 
surements were reported in Ref. by measuring the 7 decays in coincidence with 1( 0 
particles scattered inelastically from 120 Sn. These data were successfully explained by using 
the quasiparticle representation of the PDM p 6 j incorporating the thermal pairing corre¬ 
lations. These calculations are based on the modified BCS approach, where the pairing 
gap does not vanish abruptly, but decreases slowly with increasing T m EMS]. In the 
conventional BCS approach, the pairing gap collapses at a critical temperature of the tran¬ 
sition from the superfluid phase to the normal one, which is equal to T c ~ 0.57A(0), with 
A(0) being the pairing gap at T = 0. In another important low-T measurement, the GDR 
width in 1,9 Au [ 31 ] at T = 0.7 MeV was observed to be similar to the ground-state value, 
in contrast to the expectation of a larger width due to the thermal excitations. The shell 
effects were predicted [ 31 ] to explain this quenching, like in the case of 208 Pb. However, in 
Ref. [35] it was reported that a proper inclusion of shell effects leads to an increase in GDR 
width in 179 Au and there is a need for an exact treatment of the pairing fluctuations (along 
with the TSF) to overcome this discrepancy. Subsequently, the importance of considering 
pairing in the TSFM has been reported in Refs. msgsT]. 

Very recently, GDR measurements in the low-T regions were carried out at the Vari¬ 
able Energy Cyclotron Centre, Kolkata and highlighted the interesting nuclear 
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properties at low T. In a recent work HU > it has been mentioned that it would also be 
interesting to compare the data with TSFM calculations by including the effect of thermal 
pairing. In some of these recent works MU, the empirical parametrizations (which mimic 
the results of TSFM at higher T ) were reported to fail in explaining the experimental data 
at low T. It has to be noted that, at low T, the application of even the proper TSFM is 
incomplete and hence the corresponding parametrizations would also fail. This difficulty in 
such a parametrization was overcome [10] by introducing another empirical entity, i.e., the 
GDR width at critical T (which corresponds to the width observed at lowest T), and the 
GDR width below that critical T was assumed to be the same constant value ms- 

In the present work we discuss a proper TSFM that is also applicable at low T. The 
success of a modified pairing approach [TTj, [16, 26, E, 33] at low-T and that of the TSFM 
elsewhere P31I2ZIE3EBOT have motivated us to consider a combination of pairing correla¬ 
tions within the TSFM. Preliminary results of our approach can be found in Refs. PH 1361137]. 
A short discussion on our approach and some important results for the nuclei 9 'Tc, 120 Sn, 
and 208 Pb were reported recently [T8]. Here we present for the first time our complete for¬ 
malism in detail along with elaborate discussions of our results for the nuclei 9 'Tc, 120 Sn, 
and 208 Pb. Additionally, we present our results also for the nucleus 179 Au. 


II. THEORETICAL FRAMEWORK 

The theoretical approach we follow is explained in four parts corresponding to (A) the 
finite temperature Nilsson-Strutinsky method for deformation energy calculations, (B) the 
BCS approach for pairing, (C) the macroscopic approach for GDR which relates the shapes 
to GDR observables, and (D) thermal fluctuations in finite systems. 

A. Finite temperature Nilsson-Strutinsky method 

For deformation energy calculations, we adopt a formalism based on the finite tempera¬ 
ture Nilsson-Strutinsky method [T3j, which is extended to include thermal pairing. Regard¬ 
ing Strutinsky’s (microscopic-macroscopic) prescription [3T], the total free energy (Ftot) at 
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a fixed deformation is calculated by using the expression 

Ftot = -Eldm + E SF - t 1 ) 

P,N 

The sum at the right-hand side runs over protons (P) and neutrons (N). The macroscopic 
part, i.e., the liquid-drop energy (Pldm), is calculated by summing up the Coulomb and sur¬ 
face energies [32103] corresponding to a triaxially deformed shape defined by the deformation 
parameters (3 and 7 (|4j) as 

Eldm((3, 7) = {[Bs(Pn) ~ 1] +2 X [PcCS,7) - 1 }a s }A 2 '\ (2) 


The function -B s (/3, 7 ) gives the dependence of the surface energy on shape and is equal 
to the dimensionless ratio of the surface area of the shape in question to the area of the 
original sphere. The function B c ([3, 7 ) is the dimensionless ratio of the electrostatic energy of 
a distorted shape distribution to that of the sphere. Both B s (f3, 7 ) and B c ((3, 7 ) are elliptic 
integrals which are evaluated for a given deformation. The parameters a s and x are chosen 
to be 19.7 MeV and respectively, where Z and A are the charge and mass numbers of 
the nucleus. It has to be noted that E ldm(/3, 7 ) is the binding energy relative to that of 
a sphere and hence is termed the deformation energy. Such a quantity is enough for us as 
we study only the variation of energy over the deformation space rather than its absolute 
value. At finite T, in principle, one should consider that the free energy and the liquid drop 
parameters can have a T dependence. However, any effect of T on the liquid drop energy 
is found to be negligible at lower T (< 2 MeV) |44j . This will be further attenuated in the 
deformation (relative) energy that we are interested in. At higher T, the shell corrections 
vanish and hence the equilibrium shape becomes spherical where the deformation energy is 
zero. 

The calculation of the microscopic part of the free energy ([!]) is based on the Nilsson 
Hamiltonian given by (42, 36], 

Hp — + C l .1? + D( l -2(1 )) , (3) 


i 


k =1 


where the index i represents the sum over all single particles. The three oscillator frequencies 
ooki where k — 1,2, 3, are given by the Hill-Wheeler approximation m as 


uj k = uj 0 exp 



(4) 
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with the constraint of constant volume for equipotentials: oj x co y co z = Wq= constant. The 
oscillator frequency is chosen as, hu 0 = 77 MeV and the Nilsson parameters hi and /i 

values are chosen for appropriate mass regions from Refs. [351 138 , 3D] . The Hamiltonian (J3]) 
is diagonalized in cylindrical representation up to first twelve principal harmonic oscillator 
shells using the matrix elements given in [50] to obtain the singleparticle energies (e*), 
separately for protons and neutrons. 

While considering the pairing fluctuations (PF), we use the grand canonical ensemble 
(GCE), where the particle number fluctuations are allowed [131 [5Tj by fixing the chemical 
potential (A). The free energy corresponding to the thermal average within the GCE is 


F = {H 0 ) - XN P - TS , 


( 5 ) 


where Hq is the nuclear Hamiltonian (3) which is independent of T, N p is the particle 
number, and S is the entropy. The above expression can be rewritten as 

F = Y J {e i -\-E l )-2TY J Ml + exp(—Ej/T)] + ^ , (6) 

i i 

where E % = \J (e* — A) 2 + A 2 are the quasiparticle energies. A is the pairing gap, G is 


the pairing strength and they are determined using the procedures mentioned in Sec. IIB 


It has to be noted that even in self-consistent Hartree-Fock calculations ra, for a fixed 
deformation, the T dependence of e* is negligible. 

In the Strutinsky way, the smoothed free energy can be written as m 

F = 2 — A)f^ - 2T y^Sj-y^Ajkj 


+27s / f( x ) x E rii(x)dx 


( 7 ) 


x 


Here f(x) is the averaging 

( 8 ) 


with the last term included to give better plateaus conditions 
function given by 

1 q 

f(x ) = —f= exp(—x 2 ) C rn H n 

A/ 7T z ' 

m— 0 

C m = (—l) m ' /2 /[2 m (m/2)!] if m is even and C m = 0 if m is odd; x = (e — e^)/ 7 S , % is the 
smearing parameter satisfying the plateau condition dF/d^ s = 0; q is the order of smearing 
and H m (x ) are the Hermite polynomials. The averaged occupation numbers, single-particle 
entropies, and pairing numbers are given by 

poo ^ 

n t = f(x) rii(x) dx , (9) 
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Si = 


ki = 


f{x) Si(x) dx , 


( 10 ) 


/ —OO 

poo 


f(x) ki(x) dx , (11) 

J — OO 

respectively. These integrations and the one in Eq. (Jr]) are carried out numerically. In such 
a case there could be numerical uncertainties at very low T (<0.2 MeV). To avoid such 
problems it would be nice to extend the maximum term approximation method [T3] [53] 
to include pairing. The quasiparticle occupation numbers resulting from both thermal and 
pairing effects are given by 

1 


Hi = 


i — A ( Ei 

—— tanh ( — 

Ei V 2 T 


( 12 ) 


and the quasiparticle occupation numbers owing to thermal effects alone are given by 

1 


T 

n, = 


1 + exp (Ei/T) 
so that the total entropy can be written as 


(13) 


s = 2 E Si = “ 2 E K In nf + (1 - nf) ln(l - nj )] . 

i i 

The pairing numbers are given by 


(14) 


ki = -tanh (— l - 

1 2 Ei V 2 T 


(15) 


Calculating the smooth part of ^5 in this way, we come across unrealistic results such 
as a large proton pairing gap for the closed-shell nucleus 120 Sn, and a subsequent well 
deformed (/3 ~ 0.2) equilibrium shape at T = 0.1 MeV. Inherently, the Strutinsky method 
to incorporate pairing leads to an overestimation of pairing gap as discussed in Sec. Also, 


such a method leads to inconveniences while considering pairing fluctuations [Eq. (46)] where 
the calculations have to be carried out for a given pairing gap. Hence, we replace the third 
term at the right-hand side of Eq. (Jr]) with to obtain the following expression for the 
smoothed free energy: 

~ A 2 

F = 2^(ej - Xfhi - 2T^Si- — 


+ 27 s / f{x)x E rii(x)dx . 


(16) 
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In Eq. (16) the term ^7 is just a residual correction and the dominant contribution is 
through the quasiparticle energies, occupation numbers, and entropies, which are taken care 


of exactly. 

To quantify the role of pairing, we also carry out the calculations without pairing (A = 0) 
by considering the canonical ensemble (CE). In this case, the free energy is given by [T3j 

F = (H 0 ) - TS (17) 


and the corresponding discrete and smoothed energies can be written as 


F = 2Y,e,nJ -2TY, 


where 


T 

n‘ = 


1 + exp[(e* - A )/T)\ ' 

The Strutinsky smoothed free energy is given by 

F = 2 27^3, 

i i 

poo ^ 

+2q s / f(x)x nf (x)dx . 


Here we have 


~T 

n' = 


f(x) nj (x) dx 

J — OO 

and all the other quantities are the same as given in the previous section. 


(18) 


(19) 


( 20 ) 


( 21 ) 


B. BCS approach for pairing 

Here we discuss briefly how we utilize the pairing prescriptions in our formalism following 
the BCS approach At finite temperature, the superfluid phase can also exist and in some 
cases the reentrance of pairing can occur at excited states at nonzero angular momentum. 
These features have been well discussed in Refs. [5H [55j |56]. The BCS equations at finite 
temperature can be written as 

A = l-2nf) , (22) 

i 

N p = 2 - 2 n i) v ih 

i 
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(23) 




where 


v i — ~ 

2 


1 - 


(ej - A) 

Ei 


(24) 


represents the BCS occupation number and vj + uf — 1. Pairing can be dealt with either 
the constant gap approximation or the constant strength (G) approximation. In the former 
approach, A is chosen to be a constant value usually estimated from the odd-even mass 
differences 153 or taken as the average empirical value Ylj\[A MeV. In such a case only 


Eq. (22) has to be solved to obtain A at T = 0. For calculations without pairing fluctuations, 


we follow the constant-G approach with the value of G evaluated at T = 0 and solve Eqs. (22) 


and (23) simultaneously. When we consider the pairing fluctuations, we solve only Eq. (23) 


for a given A which is a variable of integration [Eq. (j46j)]. The value of G can be obtained 
from either Strutinsky calculations [58] or by fitting with empirical data [49] as outlined in 
the next section. 

Our excited state can have pairs in thermal equilibrium, which replaces the ground-state 
(T = 0) excited pairs and single particles in excited states. The pairs in thermal equilibrium 
follow the distribution defined by v?, if we neglect the influence of excited pairs and excited 
single particles. The excited pairs and excited single particles follow the distribution defined 


by nj (13). Hence we can assume that the resultant distribution is the combined effect of 
these two distributions. 


1. G from Strutinsky calculations 


In this method [59], the pairing gap corresponding to a smooth (Strutinsky smeared) 
distribution of single-particle states is assumed to be the same value obtained from the 
empirical average pairing gap, i.e., A = MeV. Then the pairing force strength can be 
calculated using the expression 


where 




( 1 / 2 ) 


+ 1 




(25) 


(26) 
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is the average density of pairs at the Fermi surface. The average level density at the Fermi 
level [g(A)j can be calculated using the relation 

1 q 

g(e) = — 7 = V'expt-x 2 ) Y] C m H m (x ) . (27) 

7 yjTT L Z ' 
i m=Q 

The Fermi energy (A) can be calculated by solving the equation n( A) = N p , and the average 
particle number (considering the states from bottom to a given energy e to be filled) can be 
evaluated using the relation 

n(e) = | \ t 1 + erf ( x )] “ 7 exp(-x 2 ) ^2 CrnHm- 1 (x) 1 . (28) 

i \ v m= 1 ) 

In these calculations, the single-particle states comprise all the states below the Fermi level 
and an equal number of states above it. 

2. G from empirical data 

In an alternate approach, the pairing strength constant can be calculated from an empir¬ 
ical fit to odd-even mass differences and their l/\/~A dependence. Such a fit has been carried 
out in Ref. |TT|] for the modified oscillator potential, yielding the values. 

Gp t N = [19.2 ± 7.4(IV - Z)]/A 2 . (29) 

In these calculations \/l5 Z and \/l5 N single-particle states above and below the Fermi level 
are included. It has to be noted that the results depend on this choice of configuration space. 


C. Macroscopic approach for GDR 


We follow a macroscopic approach to relate the nuclear shapes with the GDR observables 
[131 EH EE]. In this formalism the GDR Hamiltonian could be written as 


H = H osc + H t 


int j 


(30) 


where H osc stands for the anisotropic harmonic oscillator Hamiltonian and Hi nt characterizes 
the separable dipole-dipole interaction given by 

2 


Hint V ^ ^ 


k=x,y,z 


muJ k 

~2A 


^r^Xkiy) 


v =1 


(31) 
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where is the third projection of the Pauli isospin matrix 


A = 


1 0 
0 -1 


( 32 ) 


and q is a parameter that characterizes the isovector component of the neutron or proton 
average held 


V { N,P)(") = f 


1 ± rj- 


N-Z 


A 


u l x l( u ) ■ 

k=x,y,z 


(33) 


In terms of the dipole operator D, the GDR Hamiltonian (30) can be rewritten as 


H = //.+ 'ill'll 


(34) 


Including the pairing interactions in a simple way, the above equation can be modified as 

H = H osc + qD^D + X V ] V, (35) 

where x and P denote the strength and operator of the pairing interaction, respectively. In 
a simple oscillator description, the only relevant role of the pairing interaction is to change 
the oscillator frequencies [cu° sc {y = x,y,z)\, resulting in the new set of frequencies 

uj u = u° sc - xu V , (36) 


where 


v 

c o = 


fZAp + NA 


N 


V Z + N 1 ' ^ 

Here x having the units of MeV -1 has to be determined empirically and hence the associated 
sign is just notional with an expectation that the pairing will dampen the oscillations. Alter¬ 
natively, the role of pairing can be conceived as to renormalize the dipole-dipole interaction 
strength such that, 

V = Vo ~ XoVfuj v , (38) 

with y 0 having the units of MeV -5 / 2 . This assumption is based on the fact that the dipole- 
dipole interaction in its separable form will change only the oscillator frequencies. The 
pairing interaction also does the same and hence it is possible to combine these two effects 


as in the above relation. If pairing has to be introduced as per Eq. (36) then the parameter 
rj has to be renormalized accordingly. The u] v depends on T as inherited from A p and Apr. 
Hence when pairing vanishes, we have q = r/ 0 , allowing us to retain the parametrization 
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of 77 as in our previous approaches without pairing, and hence a meaningful comparison 
between these two approaches is possible. In the case of having nonvanishing A at larger 
T, an asymptotic value, say A ~ 0.75 MeV (T — 2 MeV) has to be subtracted from the A 
appearing in Eq. (37), to retain the asymptotic behavior ( 77| T>2 MeV = 77 o)- The presence of 


T in Eq. (38) ensures 77 | T=0 = 7/0 at any A. The exponent of T is chosen on empirical basis 
which can be improved further if more low-T data are available. 

Including the dipole-dipole and pairing interactions, the GDR frequencies in the labora¬ 
tory frame are obtained as 

oj z = (1 + 77) 1/2 o; 2 , (39) 


^2,3 = Ul + V) 




±5{(1 + <?) 2 H 2 -^) 2 } 1/2 | 


1/2 


(40) 


2 ' " y " y ' x/ J j ' 

The above relations look exactly same as the nonrotating limits of the GDR frequencies given 
in Refs. fl3, ,35; 00, 60], but here 77 is redefined to include the effect of pairing. From these 
GDR frequencies (energies), the GDR cross sections are constructed as a sum of Lorentzians 
given by 


a(E,) = ;r 


(Jr 


(41) 


i+(y-yj 2 /yr? ’ 

where E m , a m , and T are the resonance energy, peak cross section, and full width at half 
maximum, respectively [601162] . Here i represents the number of components of the GDR 


which is 3 as given by Eqs. (39) and (40). R is assumed to depend on the centroid energy 
through the relation [UP] 

R = r 0 (Ei/E 0 ) s , (42) 

where To and Eq are the resonance width and energy in the case of a spherical nucleus. 
In Ref. [63], an empirical fit between R and E t was carried out for the components of the 
ground-state GDR in several nuclei, yielding 6 — 1.9 ± 0.1 and T 0 /Eq = 0.026 ± 0.005. 
Hence, the energy dependence of GDR width can be approximated as 


R « 0.026R 


1.9 


(43) 


The value of 5 can be fine tuned to obtain a better fit with the experimental data. We 
obtained optimal fits with 6 = 1.8 for 9l Tc and 6 = 1.9 for 120 Sn, 1 ' 9 Au, and 208 Pb nuclei. 
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Nucleus 

Vo 

Xo(MeV -5 / 2 ) 

97 Tc 

2.60 

1.7 

120 Sn 

2.60 

3.5 

179 Au 

2.25 

3.0 

208pp 

3.40 

2.5 


TABLE I: The parameters r?o and xo used for different nuclei in our calculations. 


The peak cross section a m is given by p2 


2 NZ 1 . 

cr m = 60-— - (1 + a) , 


(44) 


7T a r 

where T is the full width at half maximum of the GDR cross section. It may not be 
appropriate to generalize Eq. fl44| for individual components. But it will not be possible to 
fix the three peak cross sections (cr m ) from one constraint, which is the GDR energy-weighted 


sum rule. Hence, Eq. (44) just sets the dependency of a mi on T* with a proportionality 
constant a which is given by the sum rule. The parameter a is fixed at 0.3 for all the nuclei. 
In most of the cases we normalize the peak with the experimental data and hence the choice 
of a has negligible effect on the results. The other parameters t/o (or 77) and xo (or x) vary 
with nuclei so that the experimental ground-state GDR width is reproduced. The choices 
of the parameters ? 7 0 and Xo f° r different nuclei are given in Table |T| The effects of quantal 
fluctuations other than those represented by the shell effects, such as the particle-number 
fluctuations on the pairing gap, etc., are just residual and small. We assume that while they 
are present in the measured width of the GDR built on the ground state, they are included 
in our approach by adjusting the parameters 770 and xo to reproduce that measured width 
(at T=0). Weakening of such effects with T is effectively represented through the weakening 
pairing gap. 


D. Thermal shape and pairing fluctuations 

Since the nucleus is a tiny finite system, thermal fluctuations related to the appropri¬ 
ate degrees of freedom at a finite excitation energy are large. While considering the shape 
degrees of freedom (which are dominant as far as the nuclear structure is concerned), the 
effective GDR cross sections depend on the thermal shape fluctuations and they carry in- 
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formation about the probable shape rearrangements [12] at finite excitation energy. The 
general expression for the expectation value of an observable O incorporating such thermal 
shape fluctuations has the form [3S[ ES] 


/ ? / 7 ®Hexp[-F TO T(T;A 7 )/T]0 
' t,/ T ©Hexp[-F TO T(r;/3, 7 )/r] 1 ' 

with the volume element given by T>[a] = (3 A | sin37| df3 dy. 

In PF calculations, we consider a range of pairing gap values, which are closer to the 
BCS value. The chemical potential (A) is fixed at the BCS value, thus we are allowing the 
system to have fewer bound and strong bound pairs. Hence there is a finite probability to 
have pairing even at high T also. 

The pairing held is a complex quantity since it also contains a phase and its appropriate 
metric for integration is d|A| 2 . Since the free energy depends only on A, the phase of A 
can be integrated out leading to a metric of the form A dA [6B]. By including the pairing 
fluctuations, we have 


<0> s 


Jf X, Ja„ t, V[a] exp [-F t0 t(T; ft 7 , A P , A N )/T] O 


p ,A ]\ 


In I, /a, Ja„ ®[“] ex P 1~Ftot(T ; A 7 , Ap, A n )/T] 


(46) 


with a volume element given by V[a] = /3 4 | sin37| d/3 d'j Ap A]y dAp dAjy. 


We perform the TSF calculations by evaluating numerically the integrals in Eq. (46) with 
the free energy and the observables calculated at every mesh point (the four-dimensional 
space spanned by deformations and pairing gaps), utilizing the microscopic-macroscopic 
approach outlined in Sec. |II A 


III. RESULTS AND DISCUSSION 

A. 120 Sn 

The Erst low-T GDR measurement was done in the case of 120 Sn and the low-T 
results [3D, E7j for this nucleus have been used as a benchmark for several theories. Our 
region of interest in this work is the low-T one where the earlier TSF calculations are not 
yielding satisfactory results [35] • As a first step, we include the pairing correlations in the 
TSF calculations within the BCS approach, but with free energies corresponding to the 
GCE approach [l4j. Without considering the PF, the pairing gaps averaged over the shapes 
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0.0 0.5 1.0 1.5 2.0 


T (MeV) 

FIG. 1: (Color online) Role of pairing (within the simple BCS approach) in calculations with 
thermal shape fluctuations in the case of 120 Sn. (a) Average pairing gap and (b) average quadrupole 
deformation parameter, as a function of temperature calculated without pairing in a canonical 
ensemble approach (CE) and in a grand canonical ensemble approach (GCE) using the pairing 
strengths obtained from Strutinsky calculations (G) and those quoted in Ref. [09j (G). 


could be obtained (45). The results of TSF calculations done in this way are presented in 
Fig. 0 where the pairing calculations are carried out as mentioned in Sec. |II B 1| (marked 
with the legend G) and Sec. IIB 2 (marked with the legend G). The first interesting point is 


the development of proton pairing (for Z = 50) even at low temperatures only when we use 
G. This is due to two reasons, namely (1) use of high pairing strength and (2) that pairing 
appears in deformed shapes, which can contribute when we average over the deformations. 

The former effect is mostly unphysical because, for closed-shell nuclei, no scattering in 
to the empty orbitals in the next major shell is possible. This effect is also an artifact of 


the choice of a large number of orbitals for pairing calculations (Sec. IIB 1). If we choose 
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FIG. 2: The variation of pairing gap for protons (solid line) and neutrons (dashed line), as a 
function of pairing strength in the case of 120 Sn. The prescriptions by Strutinsky calculations (G) 
and those quoted in Ref. [29] ( G ) are marked in the bottom axis. 


a restricted number of orbitals as discussed in Sec. IIB 2 then, independent of the choice of 
G , the pairing gap for a closed shell is always zero. 

The latter reason is interesting because it shows that in some cases, although the equi¬ 
librium shape is spherical and that of a closed shell, due to the thermal shape fluctuations 
the deformed shapes with pairing may contribute to the averaged values. In Fig. [jjb) the 
average quadrupole deformations ((/ 3)) are presented as a function of T. We can notice that 
the presence of proton pairing at low T, in 120 Sn, leads to a sharp increase in (f3) with its 
value raising from ~ 0.1 to ~ 0.2. The average deformation of 120 Sn at low T being ~ 0.2 is 
unrealistic, especially considering the fact that such a large deformation will yield a larger 
GDR width in contrast to the experimental trend. Also we note that, with G, the pairing 
sets in, even for the deformation /3 = 0.1. At this deformation the variation of pairing gap 
as a function of pairing strength is depicted in Fig. [2j where we have marked the values of 
the pairing strengths obtained from Strutinsky calculations ( G ) and also the prescription of 
Ref. [05] (G). The values obtained by using G are quite larger in comparison to those given 
by other prescription. If we consider the argument that pairing should lower the GDR width 
at low-T (as we know from the microscopic calculations for 120 Sn, for example as in Ref. [16]) 
to explain the experimental value, then the choice with G has to be ruled out. Hence, for 
further calculations, we choose the pairing strengths from Ref. [25]. With this prescription 
we do not have proton pairing in the deformed shapes either, and hence the averaged pairing 
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FIG. 3: (Color online) The probability distribution P(A) for gap parameter A of protons and 
neutrons in the case of 120 Sn with (5 = 0 at different temperatures. P( A) oc exp(— Ftot/T) 
where Ptot corresponds to Eq. ([ 5 ]) and we have normalized the peak of P( A). The simple BCS 
calculations yielding A = 0 translates (in the calculations with pairing fluctuations) into the 
scenario where A = 0 is the most probable value and the occurrence of A > 0 has a small but 
finite probability. 


gap for protons is zero at all temperatures as shown in Fig. [lja). The resulting (/3) with 
pairing get attenuated only at very low temperatures where the neutron pairing still exists. 
This lowering of (/3) indicates the lowering of GDR width as a consequence of the pairing 
effect. The existence of pairing at higher temperatures should decrease the width by those 
temperatures as well. Such nonvanishing or prolonged pairing gaps are predicted, e.g., in 
Refs. [Q, [32] by the calculations which include fluctuations in the pairing. 

As mentioned in the theoretical framework, the PF could be accounted for by assuming 
a grand canonical partition function where the particle-number fluctuation is allowed. The 
pairing gaps are averaged using the expression 

Ja p Ja w v [°\ ex P [-^tot(T ; p, 7, Ap, A N )/T] A* 


(A,) 


Ap,Aj> 


Ja p Ia n ex P [~ f tot(T ; P, 7 , Ap, A N )/T] 


(i = P,N), (47) 
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FIG. 4: (Color online) Variation of different probabilities (P) with respect to the order parameters 
in the nucleus 120 Sn for a fixed triaxiality parameter 7 = —180°. Solid and dashed lines represent 
the probability distributions for the gap parameter A of protons and neutrons, respectively, where 
we fix (3 = 0. Dotted lines represent the probability distributions for deformation parameter /3 where 
the pairing gaps are fixed at their BCS values (Ap = 0 and Ajv = 1.39). Different temperatures are 
assumed in different panels as in the inset. P (A or f3) cc exp(— Fpo^/T) where Ftot corresponds 
to the Eq. © and we have normalized the peak. In the top axis, the maximum value of /3 is chosen 
such that FpQT (ftmax ) — Ftot(P = 0) « F TO t(Ap = 2 . 0 ) — F T ot(Ap = 0 . 0 ). 

with T>[a] = Ap An dAp dApr- 

As mentioned in Ref. na, the gap parameter values obtained from the gap equation are 
to be understood as the most probable ones. The probability distribution for a given A can 
be calculated using P(A) a exp(— F tot /T) [where Rtot corresponds to Eq. (J5])] . These 
probabilities at various T for the case of protons and neutrons in 120 Sn are plotted in Fig. [ 3 } 
These results are quite consistent with previous results reported in Refs. [32], 51|, justifying 
our pairing calculations with the grand canonical partition function. In the fluctuation 
calculations at very low T, due to the Boltzmann factor [exp(— Ftot/T)], the most probable 
value (of either (3 or 7 or A) is favoured while we integrate over all possible values. Also, for 
the lowest value considered for T (=0.1 MeV, below which we face numerical instabilities) 
one would expect that the fluctuations should not play a role and hence the averaged gap 
has to match with the most probable gap. At higher T, the Boltzmann factor widens and 
allows a strong contribution from the values closer to the most probable ones. 
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Having stated the importance of both PF and TSF, now we proceed to have a comparison 
between the strengths of these two fluctuations. In Fig. [4], we plot the probability distribution 
of pairing gaps with a fixed shape (at (3 — 0) along with the probability distribution of 
deformation parameter f3 with the pairing gaps frozen at their BCS values. The neutron 
and proton pairing fluctuations are represented separately and these plots are shown at 
different temperatures. To enable comparison between the changes due to (3 and A, the 
upper limit of the f3 values in each panel are fixed by demanding equal energy difference 
with the change in A p values in the lower x axis. The free energy is more sensitive to the 
A p than Atv because of the nature of the proton closed shell. We can see that both the PF 
and TSF are equally important at T = 0.1 MeV. At T = 0.5 MeV, the PF dominate mainly 
because of the neutron pairing. At T = 1.0 MeV, the PF are relatively (when compared 
to previous case) less dominant. The PF and TSF again become equally important at 
T = 1.5 MeV, where the pairing field weakens and more excitation energy is available for 
the nuclei to sample higher lying states of larger deformation. For the same reason, at higher 
T (not shown here), eventually the TSF become increasingly dominant as the PF saturate. 
Importantly we infer that it is necessary to consider both PF and TSF in low-T calculations. 

Now it is interesting to see how these PF combine with the TSF. We carried out such 
calculations using free energies (|5j) from the grand partition function with averaging over 
the shapes as well as the pairing gaps for protons and neutrons. These results are presented 
in Fig. [5j where we see that the combined effect of PF and TSF (shown with the legend PF) 
is quite different from the results of TSF with pairing included through the BCS approach 
(shown with the legend BCS). It has to be noted that now when we include pairing, we do 
the calculations within the GCE ([5]), and when we neglect pairing, we resort to the CE (17). 
We can notice from Fig. [5|a) that, while having BCS pairing along with TSF we do not 
have the proton pairing gaps; but with the inclusion of PF the proton pairing develops and 
sustains even at T = 2 MeV. In addition to the averaging over A, now we have contribution 
from deformed shapes also due to the additional averaging over shapes (TSF). The single- 
particle energies change with deformation and at higher deformations we may not have a 
closed shell any longer. This may lead to a situation that the role of pairing is now enhanced 
even for a closed shell when we consider TSF in combination with PF. This enhancement is 
clearly seen when we compare the (A p) in Fig. [5|a). 

From Figs. |M and [5](c) it can be seen that both ( f3) and Y are quenched with the 
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FIG. 5: (Color online) (a) Average pairing gap, (b) average quadrupole deformation parameter, 
and (c) GDR width in the case of 120 Sn, as a function of temperature. The calculations done 
without pairing utilize a CE approach (CE), whereas the calculations with simple BCS pairing 
(BCS) and with pairing fluctuations (PF) utilize a GCE approach (GCE). The results obtained by 
using the liquid drop model (LDM) are also presented. Experimental values for 120 Sn, taken from 
Refs. [30l[67j, are shown by solid squares. For comparison, data for 119 Sb, taken from Ref. [9j, are 
also shown with open circles. 


introduction of pairing but with PF, V is more quenched than ((3). It has to be noted that 
the role of pairing on (/9) is only through the free energy surface whereas T is also affected 


by the term xV^V [Eq. (36)]. This additional term along with prolonged pairing in the 
case of PF leads to the quenching of T at higher temperatures when compared to the case 
of BCS. Thus, with the inclusion of PF, we can explain the experimental results very well. 
The calculations without pairing and with BCS could not explain the data for T at T = 1.0 
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T = 0.1 MeV T = 0.4 MeV T = 0.6MeV T = 0.8MeV 

120 Sn Without pairing 



120 Sn With pairing 


FIG. 6: (Color online) The free energy surfaces for 120 Sn at different temperatures are plotted 
without pairing (CE) and with pairing (GCE). The contour line spacing is 0.5 MeV, the most 
probable shape is marked by a solid red circle and the first two minima are shaded. 


and 1.3 MeV. In the GDR calculations with PF, for the resulting larger pairing gaps, the 
parameter x has to be readjusted. Instead, in a more effective way, we retain the same value 
for % but subtract the asymptotic (high T) value of A (typically ~ 0.75 MeV) from the 


calculated A entering Eq. (37). 

The prolonged pairing in the PF case plays no role in the free energy surfaces and hence 
(/3) obtained by using PF reaches the value corresponding to the one without pairing at the 
T where the BCS pairing gap vanishes (i.e., where it stops contributing to the free energy). 
In Fig. |5](b) we can notice that the suppression of (/3) in the BCS case is more than that 
given in the case with PF. This situation can be understood by examining the corresponding 
free energy surfaces presented in Fig. [6j These surfaces show minima at zero deformation 
while ignoring the pairing at all temperatures, and we get a nonzero (/ 3 ) due to the TSF. 
As T increases the TSF gains more strength leading to a larger (/3), even if there is no 
change in the free energies. In the present case, while T increases, the first two minima in 
the free energy surfaces span a larger area and hence compound the increase in fluctuations 
resulting in the increase of (/9). At T = 0.1 MeV, the inclusion of pairing shifts the most 
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probable deformation to a weakly deformed one from a spherical case. This leads to a small 
increase in (/3) at this point. At T = 0.4 and 0.6 MeV, we can see that the effect of pairing 
is to make the minimum crisp (or deep) which inhibits the shape fluctuations. Therefore, 
(/?) is suppressed at these temperatures, and at T = 0.8 MeV where the BCS pairing gap 
vanishes, the free energy surfaces with and without pairing look the same. If we include PF 
in addition to the role of BCS pairing, the TSF is enhanced and hence leads to a higher 
value of (/3) when compared to the TSF with BCS pairing gaps. We can also note that the 
free energies without and with pairing are calculated with CE and GCE, respectively. When 
there is no pairing, the relative energies, as depicted in Fig. [6](T = 0.8 MeV), do not clearly 
distinguish between these two ensembles. 

B. 179 Au 

Having analyzed the role of pairing and its fluctuation on GDR properties in a spherical 
nucleus, we proceed to study the deformed nucleus 179 Au. TSF calculations (without pair¬ 
ing) suggested a strong enhancement of GDR width at low T [35] owing to the deformation 
effects, whereas the experimental observation [MJ pointed otherwise. Here we explore how 
the inclusion of pairing and PF could change the previous interpretations. The calculated 
averaged pairing gaps in different cases for the nucleus 179 Au, are shown in Fig. 0a) as a 
function of T. We can see that pairing plays a strong role in this nucleus as depicted by 
the large averaged pairing gaps for both protons and neutrons. The simple BCS treatment 
suggests that the proton and neutron pairing gaps vanish at T ~ 0.4 and 0.7 MeV, respec¬ 
tively. While considering the PF along with TSF, the averaged pairing gaps continue to be 
strong even at T = 2 MeV, similar to the case of 120 Sn. 

The (/ 3) and T get quenched at low T as seen in Figs. [7](b) and[7](c), respectively. Similar 
to the case of 120 Sn, the decrease in (/ 3) with BCS calculations can be understood from the 
corresponding free energy surfaces which are shown in Fig. [8] for different temperatures and 
calculated without and with pairing. At T = 0.1 MeV, the calculations without pairing 
suggest coexisting shapes corresponding to spherical, oblate, and triaxial shapes. For such a 
situation allowing variety of shapes (shallow minimum or multiple minima in the free energy 
surface) to contribute, the TSF leads to an increase in the averaged values ((/3) and T). The 
inclusion of pairing in the deformation energies yields a deeper minimum corresponding to 


22 



0.0 0.5 1.0 1.5 2.0 

T (MeV) 


FIG. 7: (Color online) Similar to Fig. [ 5 ] but for the nucleus 179 Au. Experimental values for the 
spherical width (Fo) are taken from Ref. [33] and are shown with solid circles. The corresponding 
deformed widths (r) are estimated using the relation T = Fo(Egdr/ 14.2) 1,6 and are shown with 
solid squares. 

a triaxial shape. At higher temperatures also, this trend continues to exist. From these free 
energy surfaces, we can say that with the inclusion of pairing 179 Au can have only a narrow 
range of deformation, which means that there will be a quenching in the (/3) (which will be 
carried forward to T) when compared to calculations without pairing. At T — 0.4 MeV the 
narrowing of the minimum is strong and hence shows a strong suppression in (/3). At this 
T, the proton pairing vanishes and marks the beginning of the upward trend for (/3). In 
case of T — 0.6 and 0.8 MeV, without pairing, our calculations suggest a shallow bottom in 
the free energy surface spanning a wide range of deformations with the minimum at a large 
deformation with /3 = 0.2 and 7 = —120°. The corresponding calculations with pairing also 
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FIG. 8: (Color online) Similar to Fig. [g] but for the nucleus 179 Au. 


suggest a similar shape (with marginal preference to a deformed shape) because of weaker 
contribution from the pairing to the deformation energies. The (/3) from PF calculations are 
largerr than that of BCS calculations, for the same reason as discussed in the case of 120 Sn. 

All the effects discussed above, governing the variation of (/3), are carried forward in our 
results for T as shown in Figs. [T](b) and [T](c) . However, for T, there is an additional pairing 
effect in the GDR Hamiltonian. It is convenient to consider the following factors which can 
affect the T in our calculations [18J 

1. The pairing effects: 


(a) modification of free energy surfaces; 


(b) damping GDR frequencies through the term (35). 


2. The shell effects (in comparison to liquid drop behavior). 


The free energy surfaces of 1,9 Au and 120 Sn show the same trend while including the pairing 
effect. By including pairing, the minimum in the free energy surfaces becomes deeper and 
hence the TSF are attenuated. This results in the suppression of (/3) and T. Similarly factor 
1(b) listed above also has the same role (suppressing T) in both x ' 9 Au and 120 Sn. Factor 2 
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is expected to decrease the V in 120 Sn by favoring spherical shape, but this effect is almost 
negligible as we see that the results of LDM and CE are almost same in Fig. [5|c). However, 
we observe that this factor increases T in 179 Au. Due to this, the overall suppression of T 
(at lower T, in comparison with T = 0) is smaller in 179 Au when compared to that of 120 Sn. 

In Fig. 0c ), we have also plotted the experimental values [ 31 ], which are ambiguous. In 
Ref. [ 31 ], the (spherical) GDR width r 0 of 7 ' 9 Au at T = 0.7 MeV is quoted as 5.0 ± 0.35 
MeV with a deformation parameter /3 = 0.1 ± 0.1 and the (spherical) GDR centroid energy 
E 0 = 14.2 MeV. It may not be appropriate to compare this r 0 with the (deformed) GDR 
width (r) obtained theoretically. For a better comparison, we have estimated the equivalent 
T using the relation [34] V = Fo(Fgdr/-Fo) 1 ' 6 along with the experimentally suggested values 
for To, Eq, and f3 = 0.1. In a similar way we have estimated the width at T = 0 MeV from 
a given To (4.65 ± 0.15 MeV) [ 33 ], with /3 = —0.139 [68j. This conversion results in large 
error bars which encompass all the results corresponding to different calculations (without 
pairing, BCS, and PF). The Goldhaber-Tcllcr model calculation [@] with empirical inputs, 
yields a larger width (~5.7 MeV) compared to the experimental T 0 of 5.0 ± 0.35 MeV. All 
these facts raise the ambiguity while corroborating the results of theory and experiment. We 
also notice that, unlike the case of 120 Sn, the low-T data for 1,9 Au are not sufficient to have 
a constraint on the parameter x or to conclude whether a pairing approach is necessary 
or not. Thus, there is indeed a need for more experimental data on GDR width at low 
temperatures. 

C. 208 Pb 

The PF are observed to be important even in the spherical nucleus 120 Sn where we have 
a closed shell for protons. In order to see the role of PF in a doubly closed shell nucleus, 
we carried out the calculations for the nucleus 208 Pb and the results are shown in Fig. [ 9 } 
As expected for closed shells, the most probable values remain zero for pairing gaps for 
both protons and neutrons. With the inclusion of PF we obtain Unite (A) for both the 
protons and neutrons and they last up to T = 2 MeV. The marginal increase in (A) at 
lower T is due to reasons similar to those discussed in the case of protons in 120 Sn. The 
(/3) estimated with LDM shows a larger value when compared with other calculations which 
include the shell corrections. The PF are found to have negligible contribution to the (/9) 
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FIG. 9: (Color online) Similar to Fig. [ 5 ] but for the nucleus 208 Pb. Experimental values shown 
with solid squares correspond to those given in Refs. 

despite appreciable (Ajv) at lower T. At T — 0.5 MeV, (A n) ~ 1-0 MeV and still there 
is not appreciable effect of PF. For closed shells, the free energies are stiff with respect to 
A and hence the corresponding Boltzmann factor will be sharp, leading to a suppression of 
PF. All the effects observed in (/3) are reflected in T. Thus, in 208 Pb, the PF has no role on 
T and the shell effects strongly quench V at low T. This observation is in good accordance 
with experimental data available for T > 1.1 MeV. Measurements at lower T and with more 
accuracy can strengthen the arguments regarding the role of PF and shell effects. 
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FIG. 10: (Color online) Similar to Fig. [ 5 ] but for the nucleus 97 Tc. Experimental values shown with 
solid squares correspond to those given in Ref. El- The legend PF (GCE)* denotes the results 
obtained in the calculations by using the parameter 6 = 1.9. In all other calculations 6 = 1.8. 

D. 97 Tc 


Very recently, experimental data at low T for the nucleus 9 'Tc cn are reported and some 
of our results for this nucleus are presented in [18] . Here we discuss our results elaborately 
with the aid of free energy surfaces and strength functions. In 9 'Tc, the BCS calculations 
without fluctuations suggest that the proton and neutron pairing gaps vanish at T = 0.5 
MeV and T ~ 0.65 MeV respectively, as shown in Fig. [Td|(a). With both protons and 
neutrons in the mid shells, we obtain very large values for (A) from the calculations with 
the combination of PF and TSF. The (A) continues to be strong even at T = 2 MeV. 


The ( 18 ) and T get strongly quenched with the inclusion of pairing as shown in Figs. 10 b) 


and[l0[c). This quenching can be understood from the free energy surfaces shown in Fig. 11 
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97 Tc With pairing 

FIG. 11: (Color online) Similar to Fig. [6] but for the nucleus 9 'Tc. 

The calculations without pairing show a shallow minimum in the free energy surface, which 
allows a variety of shapes and enhances TSF. The inclusion of pairing strongly reduces the 
probable shapes as the area of the minimum gets narrower in the ((3, 7 ) plane. Subsequently, 
the TSF with pairing leads to reduced (/3) which is carried forward to T. The LDM free 
energies favor a spherical shape and hence yield a lower (f3) when compared with the results 
without pairing (which have only shell effects as additional contribution). Thus the shell 
effects tend to increase the (/ 3 ), but the pairing decreases (/?) at lower T. The role of pairing 
in the free energies is only up to T ~ 0.65 MeV, while including pairing without PF leads 
to quenching in (/?) and T. Hence, similar to the case of 120 Sn, only with the consideration 
of PF can we explain the experimental data. It is interesting to note that (/3)ldm < (/3)pf 
at low T whereas Tldm > Tpp. This is due to the fact that the role of pairing on F through 
the term 7 'P^'P is very dominant in relation to the role through (/3). 

In Fig. [!o|(c) we have shown with a short-dashed line the results from our calculation, 
where the parameter 6 is chosen to be 1.9. With such a choice, the T lowers appreciably 
leading to an apparent improvement in the fit. However, due to the uncertainties in the 
data, this distinction is not clear. For a better and detailed corroboration with experimental 
data nu, we calculated the GDR cross sections. The results for selected representative 
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FIG. 12: (Color online) The GDR strength functions 9 'Tc at different temperatures (T) are com¬ 
pared with the results obtained using the pairing fluctuations (PF) within a grand canonical en¬ 
semble approach (GCE). Experimental data are taken from Ref. [IT] and are shown by solid circles. 
The legend PF (GCE)* denotes that the calculations are with the parameter 5 = 1.9; in all other 
calculations <5 = 1.8. 


experimental data are shown in Fig. [T2j where the theoretical results with PF show a fair 
agreement with the experimental data. All the other calculations without PF lead to a 
significant shift in the centroid energy and hence poorly agree with the data. At T — 
1.4 MeV, there is an apparent increase in the experimental GDR centroid energy. This 
feature is inexplicable by our calculations but this discrepancy cannot be weighed due to 
large uncertainties in the data around the peak. However, if such a shift in centroid energy 
could be established precisely, it would shed more light on the role of pairing on the GDR 
centroid energy. The present uncertainties in the measured strength functions render the 
corresponding extracted width to be less reliable. For example, in all the cases shown in 


Fig. 12, even if we have increased the width of the theoretical cross sections, the overall 
agreement with data could be of the same quality. The calculations with 5 = 1.9 leads to 
cross sections which strongly deviate from the measured ones. Though the data for the T 
could not be used to validate the choice of 5, the data for cross sections help us to fix the 
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parameter 6 distinctly. With the unambiguously determined parameters, the TSFM with 
the inclusion of PF can explain the GDR width as well as the GDR strength functions fairly 
well at low temperatures. More precise data at low T can be more informative in this regard. 

It is worth mentioning that, as has been pointed out by Gervais et al. 1331 . the strong 
asymmetry of the GDR cross section predicted by the TSFM leads to a depletion of the 


strength in the low-energy region E 1 < 12 MeV. This feature is clearly seen in Fig. 12 


Including PF in this work improves the description of the GDR width at low T, but cannot 
resolve the deficiency of the TSFM with respect to the description of the GDR line shape, 
because this deficiency in the current framework of the TSFM is even stronger at larger 
excitation energies where the effect of PF is negligible or completely gone. 

So far, the TSFM has always followed a macroscopic approach to the GDR where the 
components of the GDR cross sections are assumed to have an asymmetric form. For 


example, in our approach it is given by Eq. (41) and by a similar equation in Refs. pjJCrE]. 
An addition of such asymmetric components with an energy-dependent width is expected to 
yield an asymmetric GDR cross section, even for spherical nuclei for which the experiments 
suggest a symmetric curve. Part of this discrepancy can be removed by choosing symmetric 
shapes for the GDR components, making use of the Breit-Wigner distributions, for example. 
However, all the parameters in our approach have to be readjusted along with revisiting the 
energy dependence of the width. Alternatively, combining a microscopic treatment of the 
GDR with the TSFM could resolve this issue. A more careful analysis in this regard is in 
progress. 


IV. SUMMARY AND CONCLUSION 

We have presented, in detail, our formalism to study giant dipole resonance (GDR) in 
nuclei at low temperature (T). This formalism is an extension of the thermal shape fluctu¬ 
ation model (TSFM) [[12], [13, 38J with the proper treatment of pairing and its fluctuations. 
In a macroscopic approach, the GDR frequencies are related to the geometry (shape) of the 
resonating system as in a general cas,e but additionally we consider the role of pairing. To 
compare with the measured GDR width, we consider the thermal fluctuations over the pos¬ 
sible degrees of freedom, namely the shape and the pairing field. This is achieved through a 
weighted average of the observable over the considered degrees of freedom. The weights are 
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Nucleus 

(i) Shell 

effects 

(ii) Pairing effects through Etot 

(iii) Pairing effects 

through xPlp 

Net 

effect 

T = 0.1 MeV 

T = 0.4 MeV 

T = 0.6 MeV 

120 Sn 

— 

t 

1 

1 

1 

1 

179 Au 

t 

1 

1 

1 

1 

fl 

2 08pb 

1 

— 

— 

— 

1 

1 

97 Tc 

t 

1 

1 

— 

1 

1 


TABLE II: Change in GDR width (T) at typically low temperatures (T < 1 MeV) due to (i) the 
shell effects, (ii) pairing effects through the free energy (Etot) at different T, and (iii) pairing 
effects through the attenuation of GDR frequencies in our macroscopic model for GDR [through 


the term xP^V (35)]. The symbols and f, represent increase and decrease in T, respectively. This 
change for the effect (i) is with respect to a liquid drop model; for effects (ii), (iii) and the net 
effect the change is with respect to the calculations without pairing but inclusive of shell effects. 
The symbol — represents no substantial change in T. 


given by the Boltzmann factor [exp(—F tot /T)] where the free energy (F TOT ) is calculated 
within the Nilsson-Strutinsky (microscopic-macroscopic) approach where the T dependence 
of the shell and pairing effects are taken care of properly. For calculations without pairing, 
we consider the canonical ensemble (CE) approach. 

While considering pairing, to calculate Etot, we consider the grand canonical ensemble 
(GCE) for which the particle number fluctuations are inherent. In a simple BCS approach, 
the paring gap vanishes at a critical T (as in a second-order phase transition) but when 
we consider the pairing fluctuations (PF), where the averaging over pairing gap is also 
included, the average pairing gap smoothly varies and remains strong even at T = 2 MeV. 
Such an extended superfluid phase is more pronounced while we treat PF in combination 
with shape fluctuations. Pairing not only contributes in modifying Etot but also plays a 
role in attenuating the GDR frequencies (and hence the width T) in our macroscopic model 


for GDR [through the term x'Pl'P (35)]. Hence, the pronounced and sustained pairing 
suggested by the calculations including PF leads to a quenching of GDR width in agreement 
with the experimental observations. This indicates that for small systems like the atomic 
nuclei, there are no sharp superfluid-normal phase transitions mm- 

For a precise match with the 120 Sn and 9 'Tc data, the consideration of PF is crucial, 
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whereas the present 1,9 Au data are not precise enough to arrive at a similar conclusion. 
In the nucleus 208 Pb, the shell effects are so dominating as to favor a spherical (closed- 
shell/unpaired) configuration and hence the role of PF is negligible. The factors affecting 
T, namely (i) shell effects, (ii) pairing effects through Ftot, and (iii) pairing effects through 
XP^V are summarized in Table [TT] for all the nuclei considered in this work. It is interesting 
to note that, although the factor (iii) is monotonic, there is a strong interplay between 
factors (i) and (ii) which can vary from nucleus to nucleus. More studies are required to pin 
down the possible correlations. 

We observe that the estimation of pairing force strength (G) in a Strutinsky way, by 
assuming an average pairing gap of 12/\/)4, leads to an overestimation as the corresponding 
results are not consistent with the observed GDR data. In the case of 120 Sn we have demon¬ 
strated this fact which opens up an idea of considering the precise low-T GDR measurements 
to provide a benchmark for pairing prescriptions. 
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